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Abstract 



^S) ' Context. In the framework of the study of the X-ray and optical emission in supernova remnants we focus on an isolated X-ray knot in the 
northern rim of the Vela SNR (Vela FilD), whose X-ray emission has been studied and discussed in Paper I. 

Aims. We aim at understanding the physical origin of the X-ray and optical emission in FilD, at understanding the role of the different physical 
^ ' processes at work, and at obtaining a key for the interpretation of future X-ray observations of SNRs. 

Methods. To this end we have pursued an accurate "forward" modeling of the interaction of the Vela SNR shock with an ISM cloud. We 

o 

oo 



perform hydrodynamic simulations and we directly compare the observables synthesized from the simulations with the data. 
Results. We explore four different model setups, choosing the values of the physical parameters on the basis of our preliminary analysis of the 
X-ray data. We synthesize X-ray emission maps and spectra filtered through the XMM-Newton EPIC-MOS instrumental response. The impact 
of a shock front at 6 million Kelvin on an elliptical cloud 30 times denser than the ambient medium describes well the shock-cloud interaction 
processes in the Vela FilD region in terms of spectral properties and morphology of the X-ray and optical emission. 
Conclusions. The bulk of the X-ray emission in the FilD knot originates in the cloud material heated by the transmitted shock front, but 
significant X-ray emission is also associated to the cloud material, which evaporates, as an effect of thermal conduction, in the intercloud 
medium. The physical origin of the FilD optical emission is associated to thermal instabilities. In the FilD knot the X-ray emission associated 
Q to the reflected shock front is negligible. 

C/2 Key words. Hydrodynamics - Shock waves ~ ISM: supernova remnants ~ ISM: clouds - ISM: kinematics and dynamics - ISM: individual 
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1. Introduction ature, chemical composition, morphology of the clouds, etc.). 

' However, the details of the heating processes and the physical 
■ The study of supernova remnants (SNRs) allows us to probe . u- u j . • »u u j • • . n . 

-' ^ ^ ' ^ scenario which determines the observed emission are still not 

the physical and chemical properties of the interstellar medium n j . j 

^ -' ^ ^ well understood. 

(ISM). Moreover, SNRs drive the mass and energy exchange in 



In this framework, shock-cloud interaction has been ex- 



the ISM, thus playing a fundamental role for its evolution. ' i ■ 

™ V • ■ c cxir. ■ u . • J u 1 tensiv el y studied both ana l ytically (e. g. IMcKee & Cowid 

The X-ray emission or SNRs is characterized by a clumpy 17-^71 ir- — : — " „ 7 r-TTrr , : — ~ ~ 

. , ,u . u u • » , A u f ,u U975, Heathcote & Brand 1983J) and numericallv (e. 

morphology that has been interpreted as a result or the in- „ , . „ „, , , r „ 1 L^^J 

, Z\ .u ui . J- u 1 J ■ g. Bedogni & Woodward 199a IStone & NormanI Il992l 

teraction between the blast wave expanding shock and in- C ." riL„„,IU^ „ „ \\,J„J , L ., If,' 'Zl 

, ,. ^ ™ ^ . , „ iKIein et alJll994 IXu & Stonelll995l and Ipraeile et alJl2005h 

homogeneities or the ambient medium. The impact of the — ; ' — ~ ~ — ; — z: \, , — ~° ; — . \ 

, , „ ., 1 J • n 1 J- • J and, recently, the combined effects of thermal conduction and 

shock front with ISM clouds influences the direction and ... ,., , . . , n , n 

1 ji^inxTr. • jj • 11 radiative cooling have been investigated for the first time by 

the speed of the SNR expansion and determines the heat- 1— — — ; , " _ _ . , , , , , ,. 

1 J o 1 .u / £ lOrlando et al. (2005). This work has shown that, when radia- 

ing and t he ionization of the clouds. Several authors (see , for . ; , . , ^ 

, rr "T^. — , 1 1I, r,r,-J 17! — . , lUoorJ J tive losses dominates, the shocked cloud is fragmented into 

example, IDecourchelle et alJll997l IBocchino et alJI2000t and , , , , , „ , ., , , , . , 

rr~ — ' ' , irrTTTTr : , , ■ r j cold and dense cloudlets, while thermal conduction determines 

IPatnaude et alJI2002l) have shown how the study of X-ray and . . ^ , . . , , , , 



optical observations of SNRs provides important information 
about the shock-cloud interaction process and about the physi- 
cal and chemical conditions of the post-shock medium (temper- 



the evaporation of the clouds and suppresses hydrodynamic 
instabilities. Despite these large theoretical efforts, however, 
no detailed comparison with the observations has yet been 
provided. In fact, although c omparisons with the observa tions 



Send offprint requests to: M. Miceli, are present, for example, in 'Hamilton & SarazinI (ll984') and 
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and do not include at the same time all the relevant physical 
effects. Recently, an observation-oriented num erical model 
has been presented by IPatnaude & FeserJ (l2005l) to describe 
the shock-cloud interaction in the southwestern region of the 
Cygnus Loop. However, no observable quantities are produced 
from this model and only a simple qualitative comparison 
between the model density plot and the overall morphology of 
the Ha emission is presented. 

iMiceli et alj (l2005i) (hereafter Paper I) analyzed an XMM- 
Newton EPIC observation of a small isolated knot (named 
FilD) in the northern rim of the Vela SNR. The analysis re- 
vealed a link between optical and X-ray emission and pro- 
vided information about density, temperature, extension along 
the line of sight, and chemical composition of the FilD cloud. In 
particular, the spectral analysis has shown that the X-ray emis- 
sion of FilD is described well by an optic ally- thin plasma with 
two thermal components at ~ 1 x 10^ K and ~ 3 x 10* K, 
respectively. The inhomogeneity of the surface brightness in 
the FilD region is determined by different volume distributions 
along the line of sight of the two components. 

Here we present a detailed hydrodynamic model of the 
shock-cloud interaction in the Vela FilD. Our approach con- 
sists in synthesizing observable quantities from the model (for 
example, maps of X-ray surface brightness and spectra), com- 
paring them directly with the XMM-Newton EPIC data and with 
optical observations of the FilD region. We also perform on the 
synthesized products the same analysis we carried on for the 
EPIC data. The aims of this work are: i) the proper interpre- 
tation of the X-ray observation (i. e. to understand the physi- 
cal origin of the two spectral components, to ascertain the evo- 
lutionary stage of the interaction, to obtain information about 
the density and temperature profiles of the emitting plasma); 
ii) to understand the role of the different physical processes 
(transmitted and reflected shocks, thermal conduction, radia- 
tive cooling) in the evolution of the system; iii) to obtain a key 
for the interpretation of future X-ray observations of middle- 
aged SNRs. 

The paper is organized as follows: in Sect. I2.ll we review 
the main results of Paper I that give us constraints for our sim- 
ulations; the model equations and setups are presented in Sect. 
I2.2l and in Sect. 12. 31 our results are presented in Sect.|3]and dis- 
cussed in Sect. |3 Finally, our conclusions are drawn in Sect. 

m 



2. Hydrodynamic modeling 

2.1. Constraints from the observations 

Fig. [l] shows the XMM-Newton EPIC MOS count rate images 
of the FilD region in the 0.3 - 0.5 keV, 0.5 - 1 keV, and 
0.3 - 2 keV energy bands. The bin size is 10.4" and the im- 
ages are adaptive ly smoothed to a signal-to-nois e ratio of 10 
and vignetting-corrected as in ' Miceh et alJ ( l2006l) . The XMM- 
Newton images show that shape and size of the FilD X-ray knot 
are similar in the three energy bands. In Paper I we found that 
the extension of the X-ray emitting knot is~(2xlxl)xl0"^ 
cm. 



The Vela FilD emits both in the X-ray and optical bands. 
Fig. |2] shows an Ha image of the FilD region, including the 
X-ray contour levels. A bright optical filament spans between 
the FilD regions with the highest X-ray surface brightness. 
The filament is almost aligned along the South-North direc- 
tion and, therefore, it is perpendicular to the plane of the inci- 
dent shock front. Optical filaments in SNRs are usually asso- 
ciated to the plasma behind the transmitted shock which trav- 
els through dense interstellar clouds, because they are typically 
aligned with the transmitted shock front planes; the peculiar 
orientation of the FilD optical filament is difficult to explain 
according to this scenario. 

As discussed in Paper I, the spectral analysis of the XMM- 
Newton EPIC data shows that the X-ray emission in the FilD 
region is described well by two thermal components in colli- 
sional ionization equilibrium (CIE) at ~ 1 x 10'' K and ~ 3x10* 
K, respectively. The emission measure (EM) of the cooler com- 
ponent is about an order of magnitude larger than that of the 
hotter component (Paper I, Fig. 8). We associated these compo- 
nents with two different phases of the clouds: the cloud core 
(which corresponds to the cooler component), surrounded by 
a hotter and less dense corona. The core post-shock density is 
ricore - 1-5-5 cm~^ and the corona post-shock density is about 
three times lower. Moreover, we calculated the distance of the 
FilD cloud from the center of the Vela shell and, by using the 
Sedov model ( Sed ovll959l) . we estimated that the shock veloc- 
ity in the interc loud medium is 5.7 - 7.8 x 10^ cm/s. Therefore, 
following McKee & Hollenbachl ( ll980l) . the shock temperature 
is Tsh = 4.6 - 8.6 X 10* K. We also derived the mean abun- 
d_ances of O (O/Oq k 1 .0), Fe {F^IFe^ = 0.39 + 0.05), and Ne 
{Ne/NeQ = 1 .7 + 0.2). We found that the inhomogeneity of the 
surface brightness in the FilD region is determined by different 
volume distributions along the line of sight of the two compo- 
nents. In particular, the surface brightness increases with the 
EM per unit area of the soft component. Finally, we estimated 
that FilD was shocked ~ 1 300 - 4500 yr ago. We will use these 
results as constraints for our modeling. 

2.2. The model equations 

We model the impact of a plane-parallel SNR shock on an iso- 
lated uniform cloud in pressure equilibrium with the ambient 
(intercloud) isothermal medium (Fig|3j- The plane-parallel ap- 
proach is justified because the FilD cloud is small (< 1 pc) 
with respect to the curvature radius of the Vela shell (~ 15 pc). 
The shock-cloud interaction is described by means of a two- 
dimensional hydrodynamic model. The model solves the time- 
dependent compressible fluid equations of mass, momentum, 
and energy conservation, including radiative losses from an op- 
tically thin thermal plasma and thermal conduction (consider- 
ing also the effects of heat flux saturation). The model equa- 
tions are: 

$+V-pv = (1) 



-^+V-pvv-HVP = (2) 
at 
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Figure 1. XMM-Newton EPIC MOS count rate images of the Vela FilD knot, observed in different energy bands. We smoothed 
the images adaptively to a signal-to-noise ratio of 10. The bin size is 10.4", corresponding to 3.9 x 10'^ cm, assuming a distance 
D = 250 pc. The color bar has a linear scale and the count rates range between 1.2 x 10 s ' bin ' and 6 x 10 s ' bin in the 
0.3 - 0.5 keV and 0.5 - 1 keV images, and between 2.4 x 10""* s"' bin"' and 1.2 x 10"^ s"' bin"', in the 0.3 - 2 keV images. 



^ + V ■ (/?£■ + P)v = -V ■ q - n.nHHT) 
ot 

with 



(3-) 2.3. The model setup 



(4) 



where v is the bulk velocity of the gas, t is the time, P the 
pressure, T the temperature, and uh are the electron and hy- 
drogen number density, respectively (we assume = «h), p = 
^mnriH is the mass density (yU = 1 .26, assuming cosmic abun- 
dances), A(T) is the radiativ e losses functio n per unit emission 
measure (Raymond & Smit hI 1 977l iMeweet al. 19 8^1), and q is 
the conductive flux, defined following lOalton & BalbusI lll993h 
asq - (l/Qsni + ^IqsatT^, where qsm and qsat are the classical 
dSpitzeJl 19621) and saturated cDalton & Balbus.1 19931) conduc- 



tive flux, respectively 



-k(T)VT . 



qsc 



-sign(Vr)50pc^ 



(5) 



with ^r) = 5.6 X 10"^r5/2 g-i cm"', ^ = 0.3 (in 
agreement with the values suggested by iGiulianil 1 19841 and 
borkowski et al. 1989 for a fully ionized cosmic gas), and 
where Cs is the isothermal sound speed. 

We s olved the set of equations by using the FLASH code 
jFrvxell et al. 2000) . upgraded to include thermal conduction 
and radiative losses fro m an optically thin thermal plasma 
(see lOrlando et al."2005V The FLASH code is based on a di- 
rectionally split PPM Riemann solver ( W oodward & Colella 



1984. and uses ada ptive mesh refinement (PARAMESH, 



MacNeice et al.*2000'). particularly appropriate to describe the 



moving steep gradients at the boundaries of the shock front. 
Our simulations have been performed on an IBM pSeries p4 
machine and on a Linux cluster of AMD opteron processors, 
for a total of ~ 9000 CPU hours. 



We show results of four simulations with different model se- 
tups. The setup of the physical parameters is dictated by the 
results of our analysis of the X-ray data. The relevant parame- 
ters are: the cloud/intercloud density contrast, the shock tem- 
perature, and the cloud dimension. Table ^ shows the initial 
physical parameters of our four setups. In all the setups the in- 
tercloud pre-shock temperature is lO'* K and the intercloud pre- 
shock density is 0.05 cm"^. In setup Sphl and Sph2 the cloud 
is spherical, whereas in setup Fill and F112 we consider an el- 
lipsoidal cloud, with the semimajor axis b^i aligned with the 
incident shock velocity. The Mach number M (i. e. the shock 
speed in units of the intercloud isothermal sound speed) and 
the density contrast (i. e. the ratio cloud/intercloud density) 
of our setups are Ms\e\ - 50, xs\e\ - 20 (for setup Sphl 
and Fill); Msiei = 57, ysj.e?. - 30 (for setup Sph2 and setup 
F112). According to lOrlando et"aD J2005> (Fig. 2), for these val- 
ues radiative cooling dominates over thermal conduction. This 
implies the local formation of thermal instabilities in structures 
with dimensions L > L,- = 1.3 x lO^T^/n {Lp ~ 10"' cm, for 
T - 2 X 10^ K an d n - 5 cm"^), which is the critical Field 
length rField"l965'). Notice, however, that the results reported 
in Orlando^ al. (2()05) were derived for an impact of a SNR 
shock with a spherical cloud with radius R = 3.09 x 10"^ cm, 
while our clouds are smaller. So we expect lower values of the 
characteristic length of temperature variations It in our simu- 
lations. Since the characteristic conductive time-scale is 



IPll 



2{y-\)K{T)T ' 



(6) 



smaller clouds imply a higher efficiency of thermal conduction, 
which may inhibit the formation of thermal instabilities. 
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Table 1. Physical parameters of the model setups. In all the setups the intercloud pre-shock temperature is 10 K and the inter- 
cloud pre-shock density is 0.05 cm"^^. The setups marked in boldface are those discussed in detail in this paper 



Model setup 


Cloud morphology 


Cloud Radius/ semiaxis 


Cloud density 


Cloud mass 


Shock temperature 






(10'** cm) 


(cm-3) 


(Mo) 


(10« K) 


Sphl 


Sphere 


Rci=2 


Kd = 1 


M,, = 0.035 


T„ = 4.6 


Sph2 


Sphere 


R,i = 3.09 


rici = 1.5 


M,, = 0.196 


Tsh = 6 


EUl 


Ellipsoid 


= Crf = 1, fed = 3.09 


rici = 1 


Mri = 0.032 


Tsi, = 4.6 


E1I2 


Ellipsoid 


ffld - Cci = 1, bci = 3.09 


«c/ = 1-5 


Mrf = 0.047 


Tsh = 6 



t v 



T - 






m 


■ ^ 






m i 






mm' 


■ „ 




N 



Figure 2. Vela FilD Ha emission (from AAO/UKST H-alpha 
survey). We have superimposed, in green, the X-ray contour 
levels at a count rate 2.4 times higher than the minimum, de- 
rived from the EPIC observation in the 0.3 - 2 ke V band (Paper 
I, Fig. 2, left panel). 

The simulations are performed in a cylindrical 2-D coor- 
dinate system (r, z), assuming axial symmetry. The computa- 
tional domain extends over 1.5 x 10'*^ cm in the r direction and 
over 2 x 10'^ cm in the z direction. We use reflection boundary 
conditions at r = and zero-gradient boundary conditions (for 
V, p, and p) elsewhere. The shock velocity is parallel to the z 
axis, which corresponds to the North direction in the observa- 
tion analyzed in Paper I. The post-shock initial conditions are 




main shock 



Figure 3. Schematic representation of the initial conditions of 
our simulations. The system consists of a planar shock-front 
which interacts with a spherical/elliptical (dashed line) isolated 
ISM cloud. 



given in the strong shock limit dZel'Dovich & Raizeilll967i) . 
The finest spatial resolution is ~ 1.95 x lO'^ cm, therefore we 
have more than 100 zones per cloud radius in setup Sphl and 
Sph2 and ^150 and ^ 50 zones per semimajor and semiminor 
axis, respectively, in setup Bill and E112. 

2.4. Synthesis of the X-ray emission 

Our modeling allows us to simulate the detailed evolution of 
the temperature and density of the shock-cloud system. From 
the computed temperature and density we are able to synthesize 
the X-ray emission detectable with the XMM-Newton EPIC- 
MOS 1 camera ( Turner et al.^2001.) . By rotating the 2D maps of 
n and T about the symmetry axis, we obtain the full 3D dis- 
tributions in the cartesian space (x', y' , z'), where the y' axis 
corresponds to the direction of the line of sight and is perpen- 
dicular to the (r, z) plane. We then derive the emission measure, 
EM{x', y', z'), and temperature, T{x', y', z'), for each fluid el- 
ement and the distribution EM{x', z') vs. T{x', z') (hereafter 
EMD) integrated along the line of sight, for each (x', z'), in 
the temperature range 0.58 - 16 x 10'' K (using 145 bins eq- 
uispaced in log T). From the EMDs, we synthesize the maps 
of the X-ray emission and the focal plane spectra using the 
MEKAL spectral code ("Mewe et al."l985', 'Mewe et al."l9M 
iLiedahl et al..l995.) . with solar abundances for all the elements, 
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but Fe (Fe/FeQ = 0.39) and Ne (Ne/Nes = 1.7), in agreement 
with the results of the data analysis (see Sect. 12. 1> . The source 
is assumed at a distance D = 250 pc, corresponding to the best 
estimate of the distance of the Vela SNR, which is known with 
a ~ 30% uncertainty (see Paper I and references therein). We 
filtered the emission through the XMM-Newton EPIC-MOSl 
instrumental response and the interstellar column density at the 
Vela FilD knot, Nh = 1 x 10^" cm"^ (see Paper I and references 
therein). The exposure time (feip'' = 50 ks, for the MOSl cam- 
era) is the same as that of the observation analyzed in Paper I 
(tf^], ~ 25 ks for each MOS camera). 

We have integrated the emission over a cube with edge 
length 10'^ cm (centered on the cloud). This length corre- 
sponds to the maximum extension of the intercloud post-shock 
medium along the line of sight and was chosen so as to match 
the minima of the X-ray synthesized surface brightness with 
those observed with XMM-Newton. 

The X-ray emission maps have been produced in the same 
energy bands analyzed in Paper I (i. e. 0.3 -0.5 keV, 0.5 - 1 keV, 
and 0.3 - 2 keV) and are presented in Sect. 13. H and Sect. 13.21 
The space binning of all the model X-ray maps is 3.9 x 10'^ 
cm that corresponds to 10.4" (i. e. the bin-size of the XMM- 
Newton EPIC images of Fig.Q at 250 pc. Focal plane spectra 
have been synthesized for different regions of the domain, to 
perform a spatially resolved spectral analysis. All spectra have 
been analyzed using XSPEC and the chosen on-axis response 
file is ml_medv9q20t5r6_all_15.rsp. 

3. Results 

3.1. The spherical cloud 

The upper panel of Fig.|3]shows, for setup Sphl, the 2D cross- 
sections through the (r, z) plane of temperature and density 
3400 yr after the first impact of the shock front with the spher- 
ical cloud. This stage of the interaction is the one which best 
approaches the observed data, as shown in Appendix A. The 
temperature and density cross-sections at different interaction 
times and the corresponding synthesized X-ray maps are avail- 
able as online material. 

As the shock interacts with the cloud, part of it is trans- 
mitted through the cloud (with temperature ~ 10^ K and den- 
sity ~ 4 cm"^) and a reflected shock propagates backwards in 
the ambient medium. Figure |3 also shows that the plasma be- 
hind the reflected bow shock has a density n ~ 0.6 cm"^ and 
it is about three times denser than the intercloud post-shock 
plasma. The pressure behind the transmitted shock front, P', 
the reflected shock front, P'', and the main shock front, f 
are: P' ^ I x 10"'' dyne cm^^^ p, ^ g 5 ^ ^q-io ^jyj^g ^^-2^ 

and P"" ~ 2.5 x 10""^ dyne cm"^. The thermal conduction be- 
tween the cloud and the intercloud medium drives a thin hot 
halo around the cloud. As the primary shock progressively en- 
gulfs the cloud, it generates lateral shocks propagating into the 
cloud. 

After ~ 3000 yr, in the regions where the transmitted shock 
interacts with the lateral shocks, there are dense (n ~ 5 cm"^), 
cold (r ~ 1 -2 X 10^ K), and elongated (L ~ 5 x 10'^ cm) struc- 
tures. These structures are larger than the critical Field length 
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Figure 4. Upper panel: temperature (on the left) and density 
(on the right) 2D cross-sections through the (r, z) plane 3400 
yr after the impact of the shock front with the spherical cloud 
(setup Sphl). Lower panel: same as upper panel for the elUp- 
soidal cloud (setup E112), 3550 yr after the shock impact. The 
color bar has a logarithmic scale. The dashed box indicates the 
field of view of the X-ray count-rate maps showed in Fig.|5l 



and are therefore thermally instable. They cool off to temper- 
atures of a few 10^ K, reaching density « ~ 10 cm"^ in a few 
10^ yr (and are visible in Fig. 0}. These values are consistent 
with the characteristic temperature and density of optical fil- 
aments in SNRs and with those observed in the FilD region 
dBocchino et al. 2000). We caution the reader that our simula- 
tions model a plasma where hydrogen is fully ionized and that 
for very low values of temperature our radiative losses func- 
tion is not calculated accurately. Therefore the values of tem- 
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Setup SpM t=i400yi 




Figure 5. Synthesized XMM-Newton EPIC MOS count rate images in different energy bands for the interaction of the Vela SNR 
shock front with a spherical cloud at f = 3400 yr after the first impact of the shock with the cloud (setup Sphl, upper panels) and 
with an ellipsoidal cloud at f = 3550 yr (setup E112, lower panels). The field of view of these image is indicated by a dashed box 
in Fig.|3] The bin size of the images is 3.9 x 10'^ cm. The color bar has a linear scale between and 1.2 x 10"^ s ' bin in the 
0.3 - 0.5 keV and 0.5 - 1 keV images, and between and 2.4 x lO"-' s"' bin"\ in the 0.3 - 2 keV images. The count rate values 
and the 2 x 10'** cm scale have been obtained assuming a distance D - 250 pc. The boxes indicates the regions for the spatially 
resolved spectral analysis (see Table|2ji. 



perature and density derived when T approaches the H ioniza- 
tion temperature may not be completely correct. However, our 
model clearly indicates that thermal instabilities can produce 
the FilD optical emission. 

The upper panels in Fig. |5l show the synthesized X-ray 
maps, for setup Sphl, at the same epoch of Fig. |5 derived in 
the three different energy bands, 0.3 - 0.5 keV, 0.5 - 1 keV, 
and 0.3 - 2 keV, by following the procedure described in Sect. 
12.41 The synthesized count rate per space bin is in agreement 
with the observed one, considering its strong dependence on 
the source distance D (the synthesized count rate scales as D^^, 
because the flux scales as and the angular extension of 
the bin as D ' x D ') and that the distance of the Vela SNR 
is known with ~ 30% uncertainty. Although the overall mor- 
phology of the synthesized emission in the 0.3 - 2 keV band 
for setup Sphl appears similar to the observed one, there are 
significant differences. Great part of the synthesized emission 
above 0.5 keV, shown in the central upper panel of Fig. |5] is 
clearly associated to the plasma heated by the reflected shock 
(see Fig. |4}, while the bulk of the soft X-ray emission in the 
0.3 - 0.5 ke V band originates in the cloud plasma heated by the 
transmitted shock. The morphology of the synthesized emis- 
sion below 0.5 keV is then very different from the synthesized 



emission in the 0.5 - 1 keV band. In setup Sphl, therefore, the 
morphology of the X-ray emitting knot in the broad band 0.3-2 
ke V is the result of the superposition of the emission associated 
with two different physical conditions of the plasma: the cloud 
material heated by the transmitted shock at North, which con- 
tributes with soft emission, and the intercloud material heated 
by the reflected shock in the central and southern parts of the 
knot, with harder emission. The extension of the X-ray emit- 
ting knot in the 0.3 - 2 keV band is significantly larger than 
that in the 0.3 - 0.5 keV band. This result does not match the 
XMM-Newton observation of the FilD knot, which shows that 
the X-ray emission in the 0.3 - 0.5 keV band has almost the 
same extension as the ones in the 0.5 - 1 keV and 0.3 - 2 keV 
bands (Fig.[l]i. 

The upper panel of Fig.|5]shows, for setup Sphl, the emis- 
sion measure of the plasma with T < 10^ K and « > 3 cm"^ 
projected in the plane of the sky, which is a proxy of the opti- 
cal emission. We also superimpose, in green, the corresponding 
X-ray contour levels at a count rate 2.4 times higher than the 
minimum, derived from the synthesized X-ray images in the 
0.3 - 2 keV band shown in Fig.|5l The impact of the Vela shock 
front with the spherical cloud of setup Sphl produces an opti- 
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Figure 6. Synthesized map (projected in the plane of the sky) of 
the emission measure of the plasma with T < 10^ K and n > 3 
cm"'', integrated along the line of sight, 3400 yr after the im- 
pact of the Vela shock front with a spherical cloud (setup Sphl, 
upper panel) and 3550 yr after the impact with an ellipsoidal 
cloud (setup E112, lower panel). As in Fig. |2] we have super- 
imposed, in green, the X-ray contour levels at a count rate 2.4 
times higher than the minimum, derived from the synthesized 
X-ray images in the 0.3 - 2 keV band shown in Fig.|5] 



cal filament aligned with the incident shock front, at odd with 
the observation (Fig.|2j- 

To study the distribution of the photon energy across the 
field of view, we synthesized the MOS mean photon energy 
map (i.e., an image where each pixel holds the mean energy of 
the detected MOS photons) in the 0.3 - 2 keV band. The map is 
shown in Fig.Q(left panel), where we superimposed the surface 
brightness contour level at 75% of the maximum in the same 
energy band. The brightest region does not correspond to the 
region with the minimum photon energy, at odd with what we 
observe in the FilD cloud (Fig.0 right panel). 

In Paper I, we singled out spatial regions with homoge- 
neous physical properties and we performed a spatially re- 
solved spectral analysis on them. In order to compare the ob- 
served spectra with those synthesized from the hydrodynamic 
simulations, we have used the same procedure, by defining 
physically homogeneous regions. The regions are indicated in 
the upper right panel of Fig.|5l for setup Sphl. In each of these 
regions there are limited fluctuations of mean photon energy 



Figure 7. Synthesized mean photon energy map for setup Sphl 
(3400 yr after the shock impact), left panel, setup E112 (3550 
yr after the shock impact), central panel, and observed mean 
photon energy map in the Vela FilD region, right panel. Each 
pixel holds the mean energy of the EPIC MOS photons in the 
0.3 - 2 keV energy band. The color scale is linear and ranges 
between 300 eV and 870 eV. We superimposed (in white) the 
corresponding contour level at 75% of the maximum count rate 
in the same energy band. 
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I T1 



channel energy (keV) 

region Y f + + A h \\ \ 
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region 7 
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Figures. Upper panel: joint fit of the observed spectrum in 
region 2 (in red) and of the synthesized spectrum for setup Sphl 
in region a (in black) with best-fit model and residuals. For the 
localization of the regions see Fig.^and Fig.|5l Lower panel: 
Same as above, for region 7 and region y. 
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Table 2. Best-fit parameters of the spectral fitting performed simultaneously on the observed MOS spectra and on those syn- 
thesized from the hydrodynamic simulations for the spectral regions indicated in Fig. [J and Fig.|5l The spectra are described 
with two thermal components (at temperature Tj and Tu) in coUisional ionization equilibrium, with an energy-independent mul- 
tiplicative factor (first row) to take into account the differences in surface brightness between the synthesized and the observed 
spectra and the different areas of the spectral regions. All the reported errors are at 90% confidence. 



Parameter (setup Sphl*) Regions a - 2 Regions a - 4 Regions yS - 4 Regions yS - 7 



0.03 



factor 0.37 ±0.02 0.45 ± 0.02 0.30 ± 0.02 0.53^ 

r, (lO^K) 0.99!0{;^ 0.99+jJ»i 1-44!J!m 1.39 ±0.07 

r„(10«K) 3.4+0.4_o.3 3.6 ±0.4 3.4+0.2_o.3 3.2 ± 0.2 

EMilEMi, 70 ±20 70 ± 20 7.2 ±1.3 4.8!°'^ 

Fe/Feo 0.4+«] 0.4 ±0.1 0.3 ±0.1 0.6!°'^ 

Reduced (d. o. f.) 1.65 (67) 1.5 (69) 4.51 (78) 2.2 (91) 

Parameter (setup E112) Regions S - 4 Regions e - 2 Regions f - 6 Regions rj -1 

factor 0.39!^-^^ 0.4 ± 0.2 l.O^jJ-^ ^■^'^Vqs 

r^lO'^K) l-ll^aos 1.04 ±0.02 ^-^""-qm 1.08 ±o!- 5 



TiiiWK) 3.5 ±0.3 3.4 ±0.2 3.3!^^ 3.5 ±0.2 

^+5 



EMilEMu 40 ±11 38+'° 24+J 8 ±3 



Pel Fee 0.5 ±0.2 0.4;0j Q.5^°^] 0.4 ±0.1 

Reduced / (d. o. f.) 1.68 (65) 1.08 (72) 2.65 (80) 2.07 (90) 

* In region y the synthesized spectrum is completely different from the observed ones (reduced > 10). 



(^ 1.7% in region a, $ 6% in region and ^ 2% in region y). 
Region a is in the bright northern part of the X-ray knot, where 
the emission associated to the transmitted shock dominates and 
the mean photon energy is low, region yS is in the brightest part 
of the cloud where both transmitted and reflected shocks con- 
tribute to the emission, and region y is located on the South, 
where we have high values of temperature and of mean pho- 
ton energy. Part of the regions we selected for the analogous 
spatially resolved spectral analysis of the XMM-Newton data 
are shown in Fig. The spectral fittings were performed si- 
multaneously on the synthesized and the observed MOS spec- 
tra. In agreement with the findings of Paper I, we adopted a 
MEKAL model of an optically-thin plasma in CIE with two 
thermal components, we fixed Nh = 1 x 10^" cm"^, and we 
left the model Fe abundance free and linked the Ne abundance 
to it, so as to have {Ne/NeQ)/{Fe/FeQ) - 4.4. We added to 
the model an energy-independent multiplicative factor to take 
into account the differences in surface brightness between the 
synthesized and the observed spectra and the different areas 
of the spectral regions. Our results are summarized in Table 
12 Notice that, since the observed spectra are by themselves 
well described by this spectral model (as shown in Paper I), the 
X'^ values in Table|2]can be considered as an indication of the 
agreement between the observed and the synthesized spectra. 
As shown in the table and in the upper panel of Fig. |S] there 
is a good agreement between the spectrum synthesized in re- 
gion a and those observed in the FilD regions with low mean 
photon energy (region 2 and region 4), but if we compare the 
spectra of region /3 and region 4 (i. e. the ones with the highest 
synthesized and observed surface brightness), we have signifi- 
cant differences (see Table |5Jl. Instead, the spectrum of region 
j6 has similar spectral features to the one observed in region 7, 
where the count rate is lower than in region 4, but the mean 
photon energy is higher. The spectrum extracted from region y, 
not shown in the table, is completely different and significantly 



harder than all the observed spectra (see the lower panel of Fig. 
|8}, even than those extracted form the RegNE cloud, which is 
the hardest X-ray emitting region in the EPIC field of view (see 
Paper I). 

In setup Sph2, the post shock temperature of the cloud (i. 
e. the temperature behind the transmitted shock) is ~ 10^ K 
(as in setup Sphl), while the cloud density is slightly higher 
than in setup Sphl (n ^ 6 cm"^). We synthesized the X-ray 
count-rate maps and focal plane spectra also for setup Sph2. 
The emission morphologies in the three energy bands (0.3-0.5 
keV, 0.5 - 1 keV, and 0.3 - 2 keV) are very similar to those of 
setup Sphl . Moreover, the synthesized X-ray emission presents 
similar spectroscopic features to setup Sphl (see Sect. 13. 1> . 
However, the global X-ray luminosity of setup Sph2 is too high 
(more than one order of magnitude) with respect to the ob- 
served one, therefore this setup will not be discussed in detail. 

3.2. The ellipsoidal cloud 

Setup Elll differs from setup Sphl for the shape of the cloud. 
However, for setup Elll, we found that in the evolution of 
the shock-cloud interaction there are no conditions for optical 
emission of the plasma'. In this case thermal instabilities do 
not develop and the plasma does not cool off to the character- 
istic temperatures of the optically emitting regions, while we 
observe a bright optical filament in the FilD cloud. Therefore, 
we will not discuss setup Elll further 

On the other hand, thermal instabilities leading to optically 
emitting regions develop in setup E112, that we discuss in de- 
tail in the following. The lower panel of Fig. 0] shows the 2D 
cross-sections through the (r, z) plane of temperature and den- 
sity ~ 3600 yr after the first impact of the shock front with the 
ellipsoidal cloud. As shown in Appendix A, this stage of the 

' We followed the evolution for ~ 6000 yr. 
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Figure 9. Upper panel: joint fit of the observed MOS spectrum 
in region 2 (in red) and of the synthesized MOS spectrum for 
setup E112 in region e (in black). Lower panel: Same as above, 
for region rj and region 7. The best fit models and residuals are 
also shown. For the localization of the regions see Fig. and 
Fig.El 



interaction is the one that best describes the observed data. The 
temperature and density cross-sections at different interaction 
times and the corresponding synthesized X-ray maps are avail- 
able as online material. 

Alt - 3550 yr, for setup E112, the plasma behind the trans- 
mitted shock has temperature T ^ 10* K and density n ~ 6 - 7 
cm^^, the plasma behind the reflected bow shock has density 
n ~ 0.4 - 0.5 cm""* and it is about two times denser than the in- 
tercloud post-shock plasma. The pressure behind the transmit- 
ted shock front, P', the reflected shock front, f, and the main 
shock front, P"", are: P' ~ 2 x 10"'' dyne cm^^, P*^ ~ 7 x 10"'° 
dyne cm"-, P'"" ~ 3 x 10"'" dyne cm"-. As shown in Fig.|4] 
the hot halo around the cloud due to its evaporation in the in- 
tercloud medium is also visible. 

The lower panels of Fig. [S] show the synthesized X-ray 
maps in the three energy bands, 0.3 - 0.5 keV, 0.5 - 1 keV, and 
0.3 - 2 keV, for setup E112. Also in this case the synthesized 
count rate per bin is consistent with the observed one (consid- 
ering a ~ 30% uncertainty for the distance of the Vela SNR), 
except for a small region in the northern edge of the synthe- 
sized X-ray knot. In the other parts of the knot the synthesized 



count rate per bin is closer to the observed one than setup Sphl . 
In setup E112 the hard X-ray emission from the plasma heated 
by the reflected shock does not contribute at large extent to the 
broadband emission. In this case, in fact, the X-ray emission 
in the 0.3 - 0.5 keV band has almost the same morphology as 
those in the 0.5 - 1 keV and 0.3 - 2 keV bands, just like in 
our XMM-Newton observation. Moreover, the overall morphol- 
ogy of the synthesized emission in all the energy bands is quite 
similar to the observed one. 

Thermal instabilities occur also in this case, in the narrow 
sheets at the lateral boundary of the cloud ~ 3000 yr after the 
shock impact. There the plasma rapidly cool off' to temperatures 
of a few lO"' K, reaching density of a few tens cm""*. The lower 
panel of Fig.|6lshows the emission measure of the plasma with 
T < IQr' K and « > 3 cm"^ projected in the plane of the sky 
for setup E112. In this case the similarity with the observations 
is more evident: the impact of the Vela shock front with the 
ellipsoidal cloud produces an optical filament aligned with the 
shock velocity, with a morphology similar to the observed one. 

Fig. (central panel) shows the synthesized MOS mean 
photon energy map in the 0.3 - 2 keV band for this setup. We 
superimposed the surface brightness contour level at 75% of 
the maximum in the same energy band. In agreement with our 
observational findings (Fig.0 right panel), the brightest region 
corresponds to the region with the minimum photon energy. 

As for the spatially resolved spectral analysis, we defined 
the four spatial regions with homogeneous physical properties 
shown in Fig. |5] (lower right panel). In these regions the fluc- 
tuations of the mean photon energy are small (^ 2% in region 
S, % 7.5% in region e and and % 5% in region 77). Region 
S covers the brightest and softest northern part of the X-ray 
knot, region e and ( correspond to the central part of the cloud, 
and region 77 is located at South, where we have higher values 
of temperature and of mean photon energy. We used the same 
spectral model described for setup Sphl to analyze simulta- 
neously the synthesized and the observed MOS spectra. Our 
results are summarized in Table |2l As shown in the table and 
in Fig. |5] there is a good agreement between synthesized and 
observed spectra. If we compare the spectra of the brightest re- 
gions (i. e. region 4 and region S) we obtain a good agreement. 
Notice that region S and region 4 are the ones with the min- 
imum mean photon energy. We obtain a satisfying agreement 
also if we compare regions with higher mean photon energy. 

4. Discussion 

We have pursued an accurate "forward" modeling of the inter- 
action of the Vela SNR shock with an ISM cloud. The direct 
comparison of the observable synthesized from the simulations 
with the data is a powerful tool to constrain the model and to 
obtain insight in the physical conditions of the Vela region. 

Our hydrodynamic model, set up from the analysis pre- 
sented in Paper I, allows us to confirm that the bulk of the X-ray 
emission in the Vela FilD knot originates in the plasma behind 
the shock transmitted into an ISM cloud. We also confirm the 
important role of the thermal conduction between the shocked 
cloud and the intercloud medium and we demonstrate that ther- 
mal instabihties can explain the peculiar optical emission ob- 
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served in the FilD region. We can ascertain the physical origin 
of the two spectral components we used in Paper I to model 
the observed X-ray spectra and we can also obtain information 
about the morphology of the ISM cloud. Our hydrodynamic 
modeling shows, in fact, how the pre-shock morphology of the 
ISM cloud influences the dynamics of the shock-cloud interac- 
tion and the corresponding X-ray and optical emission. 

Setup Sphl and E112 essentially differ for the shape of the 
ISM cloud (spherical vs. ellipsoidal), with physical parameters 
adapted so as to have similar post shock temperature behind 
the transmitted shock. Both in setup Sphl and setup E112, the 
soft X-ray emission (below 0.5 keV) is associated to the cloud 
material heated by the transmitted shock front. In setup Sph 1 a 
strong contribution to the X-ray emission above 0.5 keV comes 
from the intercloud material heated by the reflected shock and, 
in the end, the emission morphology and the spectral charac- 
teristics above 0.5 keV are not completely consistent with our 
XMM-Newton observation of FilD. In particular, as shown in 
Fig.|5j above 0.5 keV the bulk of the X-ray emission presents 
significant morphological differences from the softer emission 
(at odd with the observation) and the reflected shock emission 
is harder than the observed one (see lower panel of Fig. |8j. 
Moreover, the optical filament is parallel to the plane of the 
shock front, at odd with the observation. Fig. Qol shows the 
distribution EM vs. T of the three regions of the upper right 
panel of Fig. [S] for setup Sphl (see also Tab. |2}. In region a 
the emission measure peaks at ~ 10^ K, that is the post-shock 
temperature of the cloud material, while the contribution of the 
intercloud material at T ~ 4.6 x 10^ K is about two orders of 
magnitude lower In the central region fi the peak of the EMD 
is at ~ 1.5 X 10^ K and the contribution of the plasma with 
higher temperatures to the synthesized emission becomes con- 
siderable. In the southern region y all the emission originates in 
the hot intercloud plasma heated by the reflected shock front. 
We observe these features also for setup Sph2, which has the 
same shock temperature and cloud/intercloud density contrast 
as setup E112. Therefore, the models with a spherical cloud can- 
not explain all the spectral and morphological features we ob- 
serve in FilD. 

When the Vela blast wave shock impacts on an ellipsoidal 
cloud, the effects of the reflected shock to the X-ray emission 
are less considerable. In this case, in fact, the extension of the 
X-ray emitting knot in the 0.3 - 0.5 keV band is similar to the 
one in the 0.3 - 2 keV band (see Fig.O, just like in the Vela 
FilD cloud. Moreover, our model shows that thermal instabili- 
ties can explain the peculiar orientation of the FilD optical fil- 
ament and its relationship with the X-ray emission. As for the 
X-ray surface brightness, both models yield higher values than 
the observed ones. This result may be partially due to distance 
effects^. However, while in setup Sphl the higher count rates 
are observed above 0.5 keV, in setup E112 the brightest regions 
are those with the softest emission and the same happens in 
the FilD knot. The high surface brightness in setup E112 may 
be due to the too perfect alignment of the cloud axis and the 
main shock velocity, since the brightest knot corresponds to 

^ We remind that the distance D of the Vela SNR is known with a 
~ 30% uncertainty and that the surface brightness scales as D * 
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Figure 10. Distribution of the emission measure vs. tempera- 
ture in the spectral regions of Fig. |5] for the spherical cloud 
(setup Sphl, upper panel) and the ellipsoidal cloud (setup E112, 
lower panel) 

the region where the transmitted and the lateral shocks simul- 
taneously interact. This extreme symmetry may also explain 
why in setup E112 we observe two vertical optical filaments, 
while in the FilD region the South-North filament is only one. 
In the case of the ellipsoidal cloud, all the synthesized spectra 
are remarkably in agreement with the observation. In all the 
spectral regions the soft component dominates and the surface 
brightness is related to its EM. These results are completely 
consistent with those obtained in Paper I. The lower panel of 
Fig- El shows the distribution EM vs. T of the four spectral 
regions selected for setup E112. In all the regions the contribu- 
tion of the intercloud plasma heated by the reflected shock is 
constantly lower than that of the cloud material at lower tem- 
perature. It is worth noting that the spectral fittings performed 
on these four spectra have given results that are consistent with 
those obtained in Paper I. We conclude that setup E112 better de- 
scribes the shock-cloud interaction processes in the Vela FilD 
region. 

Our model results lead us to new physical insight for the 
interpretation of the XMM-Newton data. In Paper I we found 
that the FilD spectra are described well by two thermal com- 
ponents at T/ ~ 1 X 10'' K and T// ~ 3 x 10*' K. We associ- 
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ated these components with two different phases of the cloud: 
the cloud core and the cloud corona respectively. We can now 
better understand the physical origin of these two components. 
The cooler component clearly originates in the cloud material 
heated by the transmitted shock front. The cloud core is, in 
fact, associated to the narrow peak of the EMD shown in the 
lower panel of Fig. ^| Since the X-ray flux is always domi- 
nated by the cooler component, we conclude that in the Vela 
FilD knot the bulk of the X-ray emission originates behind the 
transmitted shock, which travels through an inhomogeneity of 
the ISM. The corona is also related to the cloud. As shown in 
Fig- El we do not have a peak of the EMD at ~ 3 x 10^ K. 
The best-fit value Tu instead reflects a wide temperature dis- 
tribution, approximatively between 1.5 x 10^ K and 5 x 10^ 
K. This is the cloud plasma that "evaporates" because of the 
thermal conduction with the hotter intercloud medium (which 
corresponds to the secondary peak of the DEM) and forms a 
hot halo around the cloud. The X-ray emission from the am- 
bient medium heated by the main shock is always negligible 
with respect to the core and corona emission. The X-ray emis- 
sion from the ambient medium heated by the reflected shock 
front is barely visible below 0.5 keV, and, although visible in 
the 0.5 - 1 keV band, is always significantly lower than the 
emission associated to the cloud. 

Although our model is focussed on a particular region of 
the Vela SNR, and we have not performed a wide exploration 
of the parameter space, we can give a few indications for the 
interpretation of future shock-cloud observations: i) the contri- 
bution to the broad band emission of the intercloud medium 
behind the main shock front is negligible in all the explored 
setups. Therefore, the X-ray emission (even the high temper- 
ature component) is associated to the inhomogeneities of the 
ISM; ii) it is possible to ascertain how much the transmitted 
and reflected shocks contribute to the X-ray emission, even if 
the statistics is not enough for a spatially resolved spectral anal- 
ysis, through the comparison between the count rate image and 
the mean photon energy map in a wide energy band: if the mean 
photon energy map (Fig.Q has a minimum in the regions with 
the highest surface brightness, the emission behind the trans- 
mitted shock front dominates, otherwise the contribution of the 
plasma behind the reflected shock is significant; iii) a typical 
approach for the determination of the plasma density of the X- 
ray emitting plasma from the spectral analysis is based on the 
assumption of pressure equilibrium between the different spec- 
tral components (this allows to derive, for each component, the 
filling factor and the density from EM and T). Our hydrody- 
namic model shows that this assumption may not be fully hold. 
However, the pressure variations are within a factor of 5, for 
the spherical cloud, and a factor of 10, for the elliptical cloud. 

We also obtained indications about the role of the different 
physical processes in the dynamics of the shock-cloud interac- 
tion. In particular, for example, the radiative losses determine 
the formation of thermal instabilities and the consequent opti- 
cal emission, while thermal conduction between the cloud and 
the intercloud medium induces the evaporation of the cloud and 
the formation of the X-ray emitting "corona". 

The efficiency of thermal conduction may be significantly 
reduced by the presence of the magnetic field (which is not 



taken into account in our model). If we assume an organized 
ambient magnetic field, the thermal conduction is anisotropic, 
because the conductive coefficient in the direction perpendic- 
ular to the field lines is several orders of magnitude lower 
than that parallel to the field lines, which coincides with the 
Spitzer's coefficient Kspi- In the presence of a pre-shock uni- 
form magnetic field, the ambient magnetic field is expected to 
envelope the shocke d cloud during the shock-cloud interaction 
(Fragile et al. 200?, where the effects of thermal conduction are 
not taken into account) and this prevents the evaporation of the 
cloud under the effect of thermal conduction and promotes the 
formation of thermal instabilities. 

However, ^ Naravan & MedvedevI 1I2OOII) have shown that, 
in the presence of a multiscale turbulent magnetic field, the 
efficiency of turbulent thermal conduction approaches the 
Spitzer's limit, the diffusion constant being ~ KspijS. Since the 
conductive flux is proportional to T^^^, this reduction of the 
efficiency of thermal conduction can be counterbalanced by 
increasing the shock temperature by a factor of 5^^^ ~ 1.58, 
and/or by slightly reducing the cloud/intercloud density con- 
trast. Since our model well reproduces the observed data, we 
conclude that there is no need to include an organized ambient 
magnetic field, but our results substantially hold in the presence 
of a turbulent magnetic field. 



5. Summary and conclusions 

Our analysis is focussed on the physical description of a bright 
isolated X-ray knot in the northern rim of the Vela SNR (Vela 
FilD), whose X-ray emission has been studied and discussed 
in Paper I. We modelled the hydrodynamics of the interaction 
between a SNR blast wave shock and an isolated cloud, taking 
into account thermal conduction and radiative cooling effects. 
We explored four different model setups, choosing the values 
of the physical parameters on the basis of our preliminary anal- 
ysis of the XMM-Newton X-ray data. Since we synthesized the 
X-ray emission filtered through the XMM-Newton EPIC-MOS 
instrumental response from the model, we can compare the 
model results directly to the data and derive detailed and quan- 
titative information and deep insight. Our study has shown that: 

- The dynamics of the shock-cloud interaction and the cor- 
responding X-ray and optical emission are significantly in- 
fluenced by the pre-shock morphology of the ISM cloud. 
In particular, if the Vela shock front impacts on a spheri- 
cal cloud, both transmitted and reflected shocks contribute 
to the X-ray emission, while in the case of an ellipsoidal 
cloud, the reflected shock contribution is negligible. 

- An eUipsoidal cloud 30 times denser than the intercloud 
medium, with the semimajor axis parallel to the shock ve- 
locity describes well the XMM-Newton observation of the 
Vela FilD knot in terms of X-ray emission spectral prop- 
erties and morphology, and in terms of spatial relationship 
between optical and X-ray emission. 

- The peculiar orientation of the FilD optical filament (which 
is difficult to explain according to the common scenario 
which associates optical emission in SNRs with slow 
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shocks travelling through dense clouds) is naturally ex- 
plained by this setup as a result of thermal instabilities. 

- Our model allows us to understand the physical origin of 
the two thermal components we used in Paper 1 to describe 
the FilD observed spectra. The cooler component originates 
in the cloud material heated by the transmitted shock front, 
while the hotter component is the result of thermal conduc- 
tion between the cloud and the hotter shocked intercloud 
medium. The heating of the cloud produces its evaporation 
and the formation of a hot diffuse halo. 

- We highlighted spectral and morphological characteristics 
of the emission from transmitted and reflected shocks, thus 
obtaining a key for the interpretation of future shock-cloud 
observations. 

Future studies should also address other open issues, like 
the presence of a single optical filament observed in FilD and 
the effects of the pre-shock symmetry of the cloud on the mor- 
phology of the optical emission and the on the surface bright- 
ness in X-rays. Moreover, we also aim at describing the shock- 
cloud interaction in RegNE, which is an other bright X-ray 
knot, hotter than FilD, discussed in Paper 1. 
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Appendix A: Determination of the evolutionary 
stage of the interaction 

The evolution of the shock-cloud interaction has been followed 
for ~ 6000 yr, for all the model setups presented in this paper. 
We have developed a procedure for the selection of the evolu- 
tionary stage of the interaction which best reproduce the ob- 
served data. This procedure is based on the definition of three 
quantities: 

1 . The extension. A, of the X-ray emitting knot, defined as the 
area covered by the region with an X-ray surface brightness 
three times larger than the minimum (in the 0.3-2 keV 
energy band). 

2. The mean photon energy, MPE, of the X-ray emitting knot, 
defined as the mean photon energy in the X-ray emitting 
region with an X-ray surface brightness three times larger 
than the minimum (in the 0.3 - 2 keV energy band). 

3. The count-rate, R, in the soft region, defined as the average 
count-rate in the region with a mean photon energy between 
the minimum £„„■„ and 1.2 x £„„„. 

Figure lA.ll shows the temporal evolution (t - corre- 
sponds to the impact of the shock front with the cloud) of A 
(upper panel), MPE (central panel), and R (lower panel), for 
setup Sphl (black diamond) and setup E112 (red cross). The 
corresponding observed values are indicated by the horizontal 
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blue lines, while the vertical green lines indicate the formation 
of thermal instabilities: no optical emission is therefore present 
on the left of the green line. 




5000 



The graphs in Fig. lA. II indicate that, for both setups, we 
can identify a temporal range, around 3500 + 500 yr where the 
synthesized values globally best approach the observed ones. 
Below f = 3000 yr we do not have optical emission, while 
above 4000 yr the synthesized values have largest deviations 
from the observed ones. For all the three quantities, the smallest 
deviations are at ~ 3500 yr, therefore in this paper we report the 
results at f = 3400 yr, for setup Sphl, and t - 3550, for setup 
E112. 

Figure IaTTI also indicates that, for t ~ 3500 yr, setup E112 
well reproduces the observed values of A and MPE, while in 
setup Sphl we have a large area of the X-ray emitting knot 
and a mean photon energy in the bright regions that is larger 
than the observed one. This is in agreement with the results 
discussed in Sect. |3 where we show that in setup Sphl the 
bright regions do not correspond to the minima in the mean 
photon energy map. As for setup E112, the high values of the 
count-rate in the soft region are due to the brightest part in the 
northern edge of the X-ray knot (see Fig.|5}, where we have the 
minimum in the mean photon energy map (Fig.0 right panel). 
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Figure A.l. Temporal evolution of the quantities A, upper 
panel, MPE, central panel, and R lower panel (defined in 
Appendix A), for setup Sphl (black diamond) and setup E112 
(red cross). The error bars are associated with the ~ 30% un- 
certainty in the distance of the Vela SNR. The blue lines cor- 
respond to the observed values and the green lines indicate the 
epoch of the formation of thermal instabilities. 



